Coherent optical control of correlation waves of spins in semiconductors 
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We calculate the dynamical fluctuation spectrum of electronic spins in a semiconductor under a 
steady-state illumination by light containing polarization squeezing correlations. Taking into ac- 
count quasi-particle lifetime and spin relaxation for this non-equilibrium situation we consider up 
to fourth order optical effects which are sensitive to the squeezing phases. We demonstrate the 
. . . possibility to control the spin fluctuations by optically modulating these phases as a function of fre- 

quency, leading to a non-Lorentzian spectrum which is very different from the thermal equilibrium 
. fluctuations in n-doped semiconductors. Specifically, in the time-domain spin-spin correlation can 

' exhibit time delays and sign flips originating from the phase modulations and correlations of polar- 

izations, respectively. For higher light intensity we expect a regime where the squeezing correlations 
(~| , will dominate the spectrum. 

=3 ■ 

1—5 ■ PACS numbers: 78.67.De, 42.50.Dv, 42.55.Sa, 42.50.Lc 

a^ ■ 

; I. INTRODUCTION 

D , When a semiconductor absorbs circularly polarized light an average collective spin polarization (S(i)) is induced in 
1-^, • the conduction band, a phenomenon known as optical orientation^^. In direct-gap /// — V semiconductors this is a 
Q ] result of optical transitions of electric dipole type across the electron-hole gap. In this process there is a net transfer of 
. angular momentum from the helicity of the photons to the angular momentum of the electrons. The optical selection 
^ ' rules are such that for the circularly polarized light the excitation rate for electrons with one spin projection is larger 
compared to the other, resulting in a net spin polarization in the conduction band. Following the post excitation spin 
dynamics one can investigate spin relaxation mechanisms using techniques such as time-resolved Faraday rotation 
and time-resolved photoluminescenceM. It is also possible to monitor the space resolved distribution of spins, their 
orientation and magnitude as well as coherently control these quantities^i^. 
^ ] Next interesting object to study are the fluctuations of the spin in a semiconductor. In thermal equilibrium these 
■ fluctuations were measured by employing the Faraday effect^. A linearly polarized probe beam that passes through the 
^ ' sample is affected by the instantaneous magnetization of the sample, and its vector of polarization acquires a rotation, 
^ , proportional to the magnetization-. By measuring the spectrum of the polarization fluctuations, it is possible to relate 
it to the spin fluctuations. 

^\ Here we propose to go beyond the equilibrium and to study dynamical fluctuations of the electronic spins which 

CS| ■ are induced by external optical held. In the way similar to the average spin induced by the polarization of absorbed 
I light we wish to consider how fluctuations of the light polarization produce upon absorption dynamical fluctuations 
\^ • of electronic spins. We will consider polarization squeezed light which has predetermined spectrum of two photon 
' correlations. The spin fluctuations spectrum will be determined by the dynamics of the absorption of such light and 
00 , will be sensitive to the phases of the optical correlation functions. These can be controlled opening possibilities for a 

' coherent control of spin correlations. 
^ , In the past theoretical suggestions to observe quantum optical effects in light-matter interaction involved models 
with squeezed radiation reservoirs interacting with atoms and semiconductoraiSiiiii^. This however requires a high 
quality squeezed reservoir, which experimentally remains challenging to achieve. Another quantum optical effect 
which does not involve a reservoir is related to two-photon absorption of squeezed light by a three-level syste m^'^i^^ 
and was demonstrated experimentally by Georgiades'^^. Other works explored the aspect of transmission of squeezed 
light through bulk media-'^^'^^'^^, and photoionization^^. Schemes for transferring correlations from light to matter 
have recently been developed in atomic and molecular optics (AMO), both theoretically^SiSi and experimentall y^^'^'^ . 
These schemes employ either coherent optical dipoles of atomic ^-systems or ground states coherence of A-systems, 
leading to a second order dependence of the spin fluctuations on the squeezed optical field. 

Unfortunately in semiconductors one is faced with strong dephasing of optical dipoles as well as valence band spins 
due to Coulomb, electron-phonon and spin-orbit interactions, rendering the above atomic optics schemes impractical. 
We note however that the dephasing and relaxation of the conduction band spins are much slower. In addition 
continuous energy bands in semiconductors have very different level structure and optical selection rules. In contrast 
to the AMO schemes, we suggest to use this and employ fourth order optical effects to manipulate the collective spin 
of the conduction band electrons through the process of two photon absorption. 

For simplicity and to isolate the transfer of fluctuations from the average we assume in this work that the semi- 
conductor is irradiated with the light which is on the average unpolarized but possesses non vanishing squeezing 
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correlations between different polarization amplitudes. Such light (polarization squeezed vacuum) was discussed by 
Karrasiov-^, Lehner— , Korolkova3^, and generated in several experiments^^'^^. In this paper we demonstrate (see Eq. 
([29| ) that by externally manipulating the phase of the frequency dependent photonic correlations a new possibility 
opens up of optical coherent control of spins in semiconductors. The essence of this effect lies in controlling the 
interference between quantum amplitudes related to many pairs of different optical frequencies. We predict that in 
order to observe this effect it is not necessary to have squeezing in quantum sense, i.e. below the shot noise limit. 
This should allow for easier measurements since such light can be generated with higher intensities. Nevertheless we 
anticipate that also in the quantum optical regime interesting features should appear in the spectrum following our 
previous predictions for the static spin correlations^^. 

The feasibility of observing spin effects in semiconductors relies on the relatively weak coupling of the conduction 
band spins to the environment, i.e. slow spin flipping processes. This has been experimentally demonstrated in 
various situations in the pasi^'^ i'^^'^^i'^^'^^ . We also rely on the fact that for the holes the situation is different with 
the corresponding rates estimated to be several orders of magnitude higher— Accordingly we have neglected the 
contribution of the hole spins to the total spin correlations. To observe induced spin fluctuations it is necessary 
to extend the experimental capabilities of measuring spin fluctuations to the non-equilibrium regime. In principle 
two beams should be employed in such an experiment, i.e. a correlated light pump and a linearly polarized probe 
acting at the same time. In addition to the measurement of the induced magnetic moment, it is necessary to have 
phase manipulation capabilities in the incoming pump beam, similarly to those demonstrated recently for squeezed 
vacuum^ ■ . 

In this work we have taken into account scattering effects of non-radiative processes in a semi-phenomenological 
way via relaxation times. However our results do not depend in an essential way on the details of the interactions 
and the main qualitative features that we demonstrate should be reproducible even in a more detailed study. This 
claim is supported by a separate diagrammatic calculation^^ which reproduces qualitatively the main result that are 
presented here. 

The structure of the paper is as follows. In ([Til introduce the two-band model of electrons interacting with the 
driving field. In sections piIIIV[) we calculate the second and fourth order contributions (in the optical field) to the 
spin fluctuations, assuming zero lattice temperature and phcnomcnological description of relaxations. In section (fV]) 
we develop physically motivated simplifications of our results, and discuss their meaning. In section (jVip we explore 
the consequences of phase modulations of the squeezing correlations on the spin fluctuations spectrum. In appendix 
Elwe give details concerning the dipole matrix elements which appear in the light-matter Hamiltonian. In appendix [Bl 
we explain the physical nature of our phenomenological approach and how it affects the results. Finally in appendix 
[C]we discuss the calculation details of additional fourth order contributions which are small and not included in the 
main calculation. 



II. MODEL 



Valence and conduction bands of a direct-gap semiconductor (/// — V) interacting with the light are modeled by 
the Hamiltonian 

^ = E ^fe^.Cfc. + '>la^k'a + Hlm + He + HgO (1) 
k{T k' a 

where the operators Cka and Vk'a- denote annihilation operators of the electrons in the conduction and valence bands, 
with quasi-momenta fc, k' and total angular momentum a. The index a in the second term enumerates both the 
degenerate valence bands and their spin degeneracy. The effect of the spin-orbit (SO) interaction in the valence bands 
is usually accounted for by Luttinger model for degenerate bands^. This model takes into account the mixing of the 
angular momentum j = 3/2 states due to spin-orbit coupling. As we explain in appendix [XI we will take into account 
the heavy-hole valence band transitions with the energies being degenerate with respect to the spin projection, and 
renormalized due to the SO interaction in the valence band. Generally important are also the Coulomb interaction 
{He) between electrons, which is responsible for excitonic and scattering effects, and spin-orbit interaction in the 
conduction band {Hso)- We will focus here on effects related to the light-matter interaction (Hlm) in the dipole 
approximation. We study optically excited electrons in the conduction band with energies above the ionization level 
of the exciton. These electrons are well described as quasi-particles with a finite lifetime arising from the momentum 
scattering induced by disorder and Coulomb interactions. We will show in appendix [X| that to a good approximation 
it is possible to write the light-matter interaction as 



k,p,a 
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where the operator Pk^k+p.a = „Ck+p,cr is the interband polarization. Q.k,p = ^k+p~^k~'^p differences between 

quasi-particle energies and the photons and 5^ ^ are optical field amplitudes with wave number q and polarization a. 
For these transitions it is sufficient to use one index (a), denoting both projection of conduction band spin (±1/2), 
projection of valence band total angular momentum (±3/2), and hclicity (±1), depending on the context where it 
appears. We also show in appendix [A] that for most purposes the interaction transition matrix element d^. ^ can be 
taken as independent of the direction of k and of s and simply denoted as a constant d. 

The spin-orbit Hamiltonian {Hso) is responsible for the D'yakonov-Perer(DP) and Bir-Aronov-Pikus (BAP) spin 
relaxation mechanisms, which are described by the following effective Hamiltonians^ 

^Scf^ = ■ a) s,s' clicks', H'^c^^^ = ^ XI X! ^)«.'''j.i'4,sCfe+g,s4',j^^fc'-gj' (3) 

s,s' k,k' ,q s,s' 

where R(k) and A are parameters depending on the material, dimensionality, temperature, and details of the optical 
excitation. We discuss in section the main effects of Hso and He on the optically excited spin correlations. 
Note that the Coulomb interaction between electrons in the conduction band (c^cc^c), which is not included, does 
not affect the relaxation of the total spin, since it is not interacting with external degrees of freedom. Even though 
in principle there are additional spin relaxation mechanisms such as Elliot- Yafet, the dominant mechanisms for 
spin relaxation are D'yakonov-Perel' and Bir-Pikus^ when the electron gas is non-degenerate (non- metallic regime). 
For observing the effects which we discuss here, materials with a relatively long spin lifetime such as GaAs are 
advantageous. Spin relaxation times of the average spin have been measured by different techniques, for different 
materials and experimental conditions^i^. These include n-GaAs quantum wells with different widths^, doping 
levels^, materials^^'^^ as well as bulk n-GaAs^. In these experiments spin life times ranging from hundreds of 
picoseconds to tens of nanoseconds have been reported. In some cases the long life times were measured even in room 
temperature?^. The contribution of the valence band holes to the driven spin fluctuations is assumed throughout to be 
small. This assumption is realistic in bulk semiconductors from the group IITV, where the spin flip rates for the holes 
are up to three orders of magnitude faster than for the conduction band electrons^, with estimates and measurements 
placing the hole spin lifetime at around Ips. These rates lead to very broad distributions of the holes spin fluctuations 
in Fourier space with negligible contributions near w = 0. It is very difficult in experiments to capture such (10'^ GHz) 
fast oscillations, so in a realistic experiment we can neglect their effect. 

Polarization properties of photons are described by the Stokes parameters^^^ which in the circular polarization 
basis are written as time averages of 

^gA(c^*)AA'VA' (4) 

, A, A' 

where cri=o..3 denote the unit and Pauli spin matrix. The averaging is over times longer than the typical correlation 
time of the field. We consider a coUinear pump beam with a range of frequencies luq ^ B/2 above the electron- 
hole gap and time and bandwidth averaged Stokes parameters (po) = {2itc/ LB)^^^Nq\, {pi) = (^2) = (pa) = 0, 
where Nq\ is the average photon occupation per mode and L is the mode quantization length. This is unpolarized 
light with the Stokes vector fluctuating around the origin of the Poincare sphere. The fluctuations are described 
by the covariance matrix pij — (piPj), which for a Gaussian type field depends on the normal {b*qyhqi\i) as well 
as anomalous {bq\bqiy) correlations, the latter constituting the main characteristics of squeezed lights. In addi- 
tion to normal correlations (&q>^&g'A') = NqdwSqq' they possess four generic anomalous correlations: two for the 

same polarization squeezing {bq±bqr±) = Mj^^Sq+qi^2qo^ujg+ujgi,2ujo and two for the opposite polarization squeezing 

{bq±bqi^) = M^'^dq+qi ^2qo5ujg+u /,2wo: whcrc m''^^ are complex functions and loq — cqq. It can be showr42iil that 

for quantum fields \M^^ \ < Nq\^Nq\ + 1 while \M^-^ '\ < Nqx for classical fields, which we discuss here. We will 
later remark on the possible effects of the quantum regime in section IVll This type of light is also called polarization- 
squeezed- vacuum2&, or two- mode squeezed vacuum state of type I^tM, 

We begin by deriving the Heisenberg equations of motion for the spin waves in the conduction band. The spin 
density operator— is given by 

5(r,t) =^V'l(r,t)a.v^.'(^,i) =EE^"'^'"''H.W^^-^'C'=^^'W (5) 

s,s' s,s' k,k' 

and we define the ^-component of the spin as 

S{q,t)^ J d\e~''^- S{r,t) = ^ cl_^.M^s,s'Ck^At) (6) 

^ k,s,s' 
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FIG. 1: (Color online) Illustration of double-photon excitation processes (long grey arrows), with yellow and red colors for 
the two pairs, which excites spin-polarized electrons (short arrows). The left (right) figure describes absorption amplitudes of 
photons of same (opposite) helicity. The interference of these amplitudes creates a spin correlation wave with a frequency Q. 



Without the optical excitation, the spin correlations are zero when the lattice temperature is T = 0, because there 
are no electrons in the conduction band. When the electric field of the light is stochastic with a zero average, so 
is the average spin component {S{q,t)), and therefore we study the correlations. As we will discuss in sections IIVI 



due to the optical selection rules only the averages {Sziq,t)Sz{q' ,t')) , with Sz{q,t) = J2k 

are affected by the squeezing correlations of the light beam (directed along the z direction) and therefore we will 
focus on them. The function {Sz{q,t)Sz{q' ,t')) is related to experiment in a similar way that the dipole fluctuations 
{a{t)^a{t')-) of an atom are related to the measurable fluorescence spectrum^"* through the flrst order correlation 
function of the optical fields^, assuming we have a suitable spectrometer at hand. In our case of the spins one way 
that comes to mind is making use of the Faraday effect, which is customarily used to measure optically injected spin 
in semiconductors. 

Using the Hamiltonian (jll2p the equations of motion which define the spin waves are given by 



dt 



''k,k' 



Ck'.a + i 



n- 



dP, 



k—p,k' ,a 



k.k' 



(7) 



dt dt 



l^^Vk',a = i^ek,k'cl^„Vk',a - id* K<y 



'k,k' -\-p,(T 



-p,k',a] + AHSO, Pl^k,a] + ^[HC, 4,fc,J 



with ^^l,k- ^4- 4'AekM = 4 - Cfc', and nl,^, 



"k.ij'- 
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We treat the interaction between the electrons and the optical field in perturbation theory. The spin operator in 
Heisenberg picture can formally be written as 

Sziq,t)^Sz{q,t)^'^ +Sziq,t)<^'^ +Sz{q,t)'^^^ +S,{q,tY'^ +Sziq,t)^^^ +... (8) 

We shall see in section Hill that to second order the spin fluctuations are proportional only to the photon occupation 
< b*b >, and therefore are not affected by squeezing, where latter only makes an effect in the fourth order. When using 
it to expand {Sz{q,t)Sz{q' ,t')) we have to take all the possible combinations which contain the optical interaction an 
even number of times, since the spin operators conserve the number of particles 

{Sziq,t)Sz{q',t')) = {Si'\q,t)Si'\q',t')) + {Si'\q,t)Si'\q' ,t')) + 

+(5(°)(g,t)5W(g',i')> + (Si''HQ,t)Si"HQ',t')) + {Si'\q,t)Si'\q' ,t')) + 

+ {Si'\q,t)Si'\q',t')) + {SP{q,t)SP{q',t')) 

(9) 

However, at T = the combinations 0x0, 4x0, 0x4 are zero since the spin operator (on the "0" side) would be 
acting on the empty conduction band state. Therefore we are left with 

{Sz{q, t)Sz (g', t')) = (5(1) (<Z, O^i'^ (g', t')) + (5(2) (q, t)5p) (g', t')) + (10) 

+ (5(i)(g,t)5f )(q',t')> + (5f )(g,t)5(i)(g',<')) 
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i.e. only 2nd and 4th order contributions. The former are not affected by squeezing and contribute to the background 
of the fluctuations spectrum. The latter, however, is affected by the q-dependent phase modulations of the function 

M^^'^\ through a microscopic process of double absorption (see section llV|) . This process is illustrated in Fig. [1] 
showing two cases of opposite and same polarization. Two quantum amplitudes related to absorption of two different 
pairs of photons proportional to M^-^'^'>{uj) and M'^^'^\ui + fi) will be shown to appear in interference terms like 
M*{w)M{uj + which bring into effect their externally controlled phase difference. 

We shall derive the Heisenberg equations of motion up to 4th order in the optical field, by writing formal solutions 
and substituting back to the previous order. Then we will use the resulting expansion to calculate the averages in (jlOp 
order by order, taking into account the nature of the optical field. An important part of the model is the treatment 
of electron and spin lifetimes, which we will discuss in section pff[) . 



III. CALCULATION OF THE SECOND ORDER CONTRIBUTION 



This physical optical process, of the lowest order, is responsible for generating carrier and spin densities in the 
conduction band of a semiconductor. It is convenient to analyze their space-time profile in Fourier space (q, lo), given 
by the correlation function of Sz{q,i^) operators. The perturbative correction which we discuss here turns out not 
to depend on squeezing or coherent properties of the light field, and will ultimately just serve as a background for 
the more interesting higher order processes (2 x 2). At the last stage of the calculation (see Eq. p9|)) we neglect 
the momenta of the spin wave {q) and the photon (p) compared with the typical electronic momentum (fc), in places 
where they appear together as a sum (e.g. k + p k). Within this approximation the only angular dependence that 
remains is that of the the dipole matrix element (see appendix [X| , which amounts to renormalization of the value 
of the matrix element. Now we use ([7]) to construct the equation for the spin, with 6^^ (bq.i) denoting right (left) 
polarizations, respectively. The first equation then reads 



dt 



E 

k 



-dPl 



''k-q,k,l/ 



E 

k 



-dP, 



k—p,k- 



i[Hso,S,{q,t)]+t[Hc,S,{q,t)] 



d*b;^Pk. 



p-q,k,] 



(11) 



The first term on the r.h.s describes a free evolution of the spin wave with typical frequencies of v^q, where vq is the 
electron velocity at fcg, the latter being the quasi- momentum of the electron optically excited by the photon of the 
middle frequency wq- Since the light induces spin correlations with the dispersion of cq, we are physically motivated 
to neglect the slow free evolution {v^q) in the following treatment. 

Instead of developing a microscopic theory for the relaxation of the spin wave Sz and interband polarization 
P, we replace the combined effect of Hso and Hq by a quasi-particle description with a width 7 = and a 
phenomenological relaxation 7^ = t~^. Under broadband excitation conditions, and when the electrons kinetic energy 
is large compared to the electron-hole exchange energy, it is reasonable to model them as quasi-particles with a finite 
life time (t). This lifetime enters the model through the level width of the electronic energies e^. In addition to r it is 
physically plausible, and supported my many experiments, to assume the existence of another, macroscopic relaxation 
time [ts) of the spin wave which is usually much longer than the electron lifetime, Ts ^ t for semiconductors such as 
GaAs. The spin lifetime can be entered as a phenomenological decay term (—7^5'^) in the equation of motion for the 
spin operator Sz (Eq. (jlip ). similarly to a decay term in a quantum-Langevin equation of motion. Equivalently, it 
can enter as an additional imaginary part (broadening) to the frequency ujq of the spin wave. We explain in appendix 
|B]the limitations of the phenomenological approach, and how to reconcile it with the results that we obtained for the 
static correlations in the previous work^^. The averaging of {Sz{q,t)Sz{q' ,t')) over the ground state of with T = 
leaves only the terms (PP^). The general expression reads 



{Sziq,t)Sziq',t')) 



(12) 



= |rf|2g-7.(t+0 rdt^e-^st, r dt.e^^'^Yl E (-l)'^~'^'e'^"''*^~"'''*^'(^;aV.')(^;c-p-,,...(ii)P2'-p',fe'-,',.'(^2)) 

■''o k,k' p,p' ,a,(7' 

where to is the initial time when the system was in the ground state before the optical fields were turned on. Given 
that Npa- is the photon occupation function, 7 is the quasi-particle lifetime, and cr,a' = ±1/2, we also have 



^p.p' f^cr,(j' A^pcr 



{bpabp'a' 



{Pk-p-q,kM{h)Py -p' . 



(13) 



-9 ,o" 



'(^2)) — S„,a'Sk-p-q,k'~p'Sk,k'-q'e 



-t2)g-7|tl 



-t2| 
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Working with this expression it is straightforward to get for t' > t, and steady state t — to ^ -f ^, 



g-(infc,p,,+7)(t'-t) 
Is - il + i^k,p,q)'^ 



-7=(*'-*) 



7s (7s 



(14) 



where Np — 'Ylia-^p,'^ ^fc,p,g — ^ + ^k-p-q- Note that since d is the interaction matrix element of the 
Ught-matter interaction, it scales like V~^^^ because of the electric field normalization. Therefore the correlations 
scale with the volume as V and are dimensionless (with our definition of Sz)- It is instructive to study the structure 
of the expression inside the brackets in Eq. p4p . First we note that it is stationary since the electronic correlations in 
Eq. (|13p and the optical fields are also stationary. For a fixed k and p we expect naturally for the frequency ^k,p,q to 
appear with t — t' (the first term). However we also expect a buildup of a constant carrier density in the conduction 
band, which carries an inevitable static spin density fluctuation (Sf). Therefore we expect the spin fluctuation to 
have a DC Fourier component (the second term), balanced by the spin relaxation 7^, and this indeed will turn out to 
be the case when we evaluate the /c-summation (see below section |V| . 



IV. EFFECTS OF SQUEEZING 

Let us now turn to the fourth order contribution to the correlations. Using expression (|Cip . we first average 
over the electronic T — ground state of the semiconductor (denoted by (•)ez)- This averaging can be done on 
Sz{q,tY^'> independent of Sz{q' ,t')'-^^ since for T = there are no connected parts between the operators n-Jl'^/ in 
{Sz{q,tY'^lSz{q',t')^'^^). This gives, for steady state (i - to > 7"\77^) 

{Sz{q,tY'^)ei = [5fc,p..=T -'5fc.p,.=i] (15) 

k,p 

with 

-i{ujp-ujp-g)tu* u i(wp-Wp+,)tL* L 
^ _ "p^qa-"p(^ I "pg-"p+i?g- ^-j^g^ 

'^''^ (7s - "i-i^p - ^p-q)){l - i{^ik-q,k-p - ^^p-q)) [is + j(Wp " UJp+q)){'^ + i{Aek,k-p-q ^ ^p+q)) 

where we introduce 7, 7^ already at the operator level as in the Heisenberg-Langevin approach, this being mathemat- 
ically equivalent to adding them at the end as representing a single particle level broadening and spin wave frequency 
broadening. For the spin correlations we have the following dependence on fourth order correlations of the photon 
fields ((•) denotes full averaging) 

{Sz{q,t)Sziq',t'))(''^ = \d\^Yl E [ (17) 

k,p k' ,p' ,X 

(7s - i{uJp - LL!p-q)){'J - i{Aek-q.k-p - LJp-q)){'Js - i{i^p' - LUp'-q')){-f - i{Aek'-q',k'-p' - UJp'-q')) 

e-("p---.+.)*e»K'-.'+,')*' [{b;,bp+q,b;,,bp,+q,,) {b;,bp+q,b;,-^bp,+q,-,)] ^ 

(7s + i{ujp - UJp+q)){j + i{Aek,k-p~q - t^p+9))(7s + ii^^p' - ^p'+q')){l + i{Aek'M'-p'~q' - LOp'+q')) 

e-^('-.-'-.-.)te^('^^'-'^^'+^'^''[{b;_^^bp^b;,^^^^^^ 



(7s - i{ujp - Wp_g))(7 - i{Aek-q.,k~p - LOp-q)){ls + i(Wp' - Wp'+q'))(7 + i{A€k'.k'-p'-q' - ^p'+q')) 

g.(.p~.p+,)*g-»(.p,-.p,_,,)*' [(6;,&p+,,fo;,_^„ Va) - %xhp+qxbi,_^,-^bp,-^)] 



+ 



(7s + iiujp - Ulp+q))^^ + i{Atk^k-p-q ~ OJp+q)){ns " «(Wp' - LOp' ^q:)){-f - i{Aek> -q' ,k' -p' " ^^p'~q')) 

where each photon correlator can be factorized into products using intensity (6*^ Ai^<2i-Ai) ^^d squeezing {bq^^x^bq^^x^) 
since the field distribution is Gaussian for light which is down-converted using a x^^^ non-linearity42. As in the case 
of static spin correlations^^, we assume a symmetric distribution of the correlations as functions of wave vectors, 
i.e. q + = Spo, which holds for non-degenerate (broad spectrum) down-converted light "^P. We can write for the 
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FIG. 2: A diagrammatic representation of the resonant interaction processes in the normal and anomalous channels leading 
to spin correlations. The normal channel (a) is related to the auto-correlation function A'^ of optical modes with themselves, 
whereas the anomalous channel (b) involve a correlation between different optical modes, either of same polarization (A/i) or 
of opposite polarization (M2). 



correlations 

(&;-g,T^P-T&;%T V-9',t) = Nj,^qNp5q^q,5p^p, + M^i);^/^!)^,,,, V+p-«,2po (18) 

(^P,T^P+'7:T^P',tV-9',T> = NpNp+gSq^q>Sp+qy + M^'^^* M^^^^Sq^q' Sp' +p,2pQ 

(^p,T^P+<7:T^p'+g',TV,T> = NpNp+q5q^q.5p..p + M^^^*M'^^^5q^q>5p>+p+q^2po 



where po is the central wave- vector of the optical spectrum, and note we omitted the ± index assuming Af^" = 
for simplicity. In figure ^ we can see a diagrammatic representation of the microscopic processes that lead to the 'N' 
and 'M' terms in the expression p8| . These are essentially two particle-hole excitations which are correlated through 
the existence of photonic correlations. These photon correlations are spectral functions, which can be externally 
controlled at the sourceSSi^liS^. Given these functions it is possible to characterize the fluctuations of the Stokes 
vector, whose average is assumed to be zero in this calculation (unpolarized)^. 

We note that there is also a (^^(g, i)*^^^S'z(g', t')*-^-*) contribution which we discuss in detail in appendix [Cl This 
contribution has a much smaller phase space compared to {Sz{q, t)'-^^Sz{q' , t')'^^) and therefore we neglect it. Also, we 
have shown^ that a direct diagrammatic evaluation of the 4th order diagrams. Fig. [51 done in the Keldysh two-time 
formalism gives a result similar to expression p7p . with the only difference being the replacement 7^ 7. The 
reason for this is that in the microscopic calculation we have not taken into account explicitly the spin relaxation 
mechanisms, which is a topic for further research. 



V. DISCUSSION AND RESULTS 



We would like now to discuss and simplify the expressions for the spin-spin correlations obtained above. We begin 
from the second order contribution (|14p which is quite a complicated integral when the vector nature of k,p,q is 
taken into account. First note that for the problem we have in mind p k, since typically p ~ lfim~^ whereas 
fco ^ for electrons with kinetic energy of iOmeV. In Fig. ^ results of numerical evaluation are shown for 

the time Fourier transform of Eq. (|14p . assuming a unidirectional beam p\\z and spin wave vectors in the direction of 
the beam q\ \p. We see that the q-dependence is very slow for q <^ ko, and similarly it can be shown that the integrand 
dependence on p is negligible for p <ti kg. Therefore a reasonable zeroth approximation would be to neglect p and q 
altogether w.r.t. k. If we further restrict ourselves to small spin waves momenta q < B/c (B is the optical bandwidth 
in units of frequency), then q <^ p and certainly this approximation is valid. We will also assume a constant electronic 
density of states in the energy range = ±7 where the integrand is appreciable and get a simple expression for the 
spin-spin correlations {per unit volume, and restoring missing ?i's) 

{Sz{q,t)Sz{q\t'))^^^ ^ 5,_,,!^^!^i^e--l*-*'l (19) 

where Pei{£o) = V'2m'^e/h'^\e=eo is the electron-hole density of states per unit volume in the bulk, taken at the central 
energy of excited electrons eo = hujo — Eg, with m the reduced e-h efl'ective mass, and e the excess kinetic energy of 
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FIG. 3: The spin fluctuations spectrum (second order contribution) plotted as a function of {q,u}) for excitations with a wave 
vector parallel to a unidirectional beam, scaled with the central electronic wave vector fco and the spin relaxation rate 7s, 
respectively. 



the e-h pair above the gap. Denoting {{S^q^w) — F.T.[{Sz{q,t)Sz{q' ,t'))] we obtain 

u^2, ,(2) 2n\d\'Peij:,N, 2n\d\^p,,Jdu;'po,,{u;')N{u;') 

\{Jz)q,i^/ — 1/ 9 , ^ — 1-/ o , ^s, 



where N{uj') = Nq\^r^cq is the photon distribution function in frequency space, and popt{i^) = w^/tt^c^ is the free space 
optical density of states per unit volume^. The contribution for long wavelengths is thus a Lorentzian background 
whose width is determined by the spin relaxation time. 

Turning now to the fourth order contribution, we observe that in certain physical situations it (|17p can be simplified 
considerably. This happens if the excess energy of the light above the electron-hole gap {huo — Eg) is large compared 
to the exciton binding energy. In this regime, using the correlations p8)) . and assuming that k ^ p, q we get for the 
only important contribution the following approximate simple expression 

{Sziq,t)Sz{q\t')y^^ = Sq^- /''';^f' Et (21) 

^ i^q + Is) p 
Np^q^ + Np,+ + Np^ + Np+q^+ + 7Vp_,^_iVp,_ + Np^^Np+q^^ + 



where pei is the electron-hole joint density of states at Ae^ = TiUi^. Since the field is unpolarized = N^, and 
assuming again for simplicity the symmetric cas 
Therefore we can omit the polarization indices 



assuming again for simplicity the symmetric case M^'X — M^^^ = Mq"'^ , to which we always refer from this point on. 



4_2 2 |J|4 -iwg(t-t') 

{Sz{q.t)Sz{q'.t'))^''^ ^Sq^^q \ , Y.^Np^^^p + NpNp+q+ (22) 

p—q p ' p p+q ^*^^p—q^^^p ^^^p p+q\ 

The restriction q' = —q which reflects translational invariance leads to {Sz{q,t)Sz{q' ,t')) — {Sz{q,t')Sz{q' ,t))* . The 
r.h.s. of Eq. (|22|) indeed obeys this if the summation over p is real, which can be shown explicitly by using the 
symmetry property Mp = M2pg-p of the squeezing correlation. 

It is useful at this point to compare the decay and oscillatory time dependence of Eq. (fTO|) and Eq. (PT|) respectively. 
The second order contribution describes part of the spin fluctuations which is only due to the effect of the 
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FIG. 4: Diagram of the second order contribution to the spin-spin correlations. The electron in the valence band (dashed line) 
is excited into the conduction band (solid line) by the photon correlation (dotted line). The only possible dynamics between 
the times t,t' in this order of perturbation theory is due to the phenomenological spin relaxation (ts). 



environment. The spin remains temporally correlated over a time scale Ts that it takes for the relaxation processes to 
be effective. As can be seen from Fig. 21 the two orders of the interaction with the light are not sufficient to affect the 
dynamics of the correlation function {Sz{q,t)Sz{q' ,t')) . The contribution of the fourth order processes is physically 
different. Here the spin correlation dynamics reflects directly the correlations of the polarization of the light. This can 
be understood also diagrammatically from Fig. [2l In the fourth order contribution there is no electronic propagator 
temporally connecting Sz{q,t) and Sz{q',t') (see Fig. [5]). Only the correlation functions of the light are connecting 
the two electronic diagrams, and the dynamics of the light is affecting the two times dynamics. We see that the 
spin relaxation 7s does affect the weight of each spectral component Sz{q,t). Fast spectral components ujq ^ 7s are 
independent of 7s and their weight is oc w^^, i.e. fast oscillations do not have enough time to be affected by the 
spin relaxation. In contrast, slow spectral components coq 7s are building up slow enough to be balanced by 7s, 
and their weight is approximately cx 77^. In particular, without the spin relaxation in the model we would have a 
divergence of the q = component, which corresponds to an unbalanced accumulation of spin-correlated charges in 
the conduction band due to the double absorption process. Note that there are also other diagrams of the fourth 
order where 7s can affect the time dependence, however they have a much smaller contribution due to electronic phase 
space considerations (see discussion in appendix [C] and elsewhere^) . 
The spectral function looks like 



^,^_ A7r^pl\d\^S{u;-cq) 



E 



(23) 



p-q p 



p+q p~q p p p+q 



The (5- function reflects the dispersion of the light. The p summation can be transformed to a frequency integral (so 
now the l.h.s. is per unit volume) 



iiS'z 



),, = 4^ ^e'M '^('^ [ dc.'popt{Lu') [N{u;' - u;)Niu;') + 7V(^')A^(^' + ^)+ 



(24) 



We see that as the phases of M^^'^) ^re modulated as a function of frequency, the spectrum of the driven spin 
fluctuations is also modulated. The spin structure factor in these situations is in fact coherently controlled using the 
phases of the w-dependent squeezing. Combining the second (|20p and fourth (|24p order contributions, defining 



C(w) = y dJ {N{J - uj)N{J) + N(uj')N(uj' + w)-h 

M(i)* (cj' - w) Af {J) -f m(i)* (cj')M'^^ (^' + ^) - 
-M(2)*(c^' - cj)A/(2)(c^') - A/(2)*(c.,')Af(2)(^' +,^) 



(25) 



and taking the optical density of states as approximately constant, given by popt — Popt{LOo) from the middle of the 
spectrum, we obtain 



/(q2\ \(2+4) _ '2n\d\'^ Pel Popt 



2Pel\d\^ 



dw'N{uj') + 4^^^l2i^5(^ „ cq)C{oj) 
n 



(26) 



Note that this result is valid for long wavelengths {q B /c) , which is also the range where the spin structure factor 
is largest. 
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Let us now choose as direction of light propagation z and integrate over qz, setting qx = % = 0. Since our 
approximations are only valid for small we will use a cutoff B /c on the q^ integration, which mean that the result is 
the w-spectrum spatially averaged over the coherence length c/ B in z-direction and the x — y plane. Physically, this is 
motivated by having in mind a Faraday probe beam passing through the sample and measuring the spin fluctuations. 
For small frequencies uj <^ B this integration will smear-out the J-function singularity. We assume that the sample 
is much smaller compared to the coherence length of the light Ic, so the averaging is effectively done over the whole 
sample. This integration gives the following spin spectral density 



uq2-. \(2+4) _ f^'" dn ((^'^\ \(2+4)| _ 47r|dpPeiPopt^ 



duj N(uj j H — C{uj) 

hB 



(27) 



In order to estimate the strength of the coherent effects compared to the background of second order fluctuations 
we need to calculate the ratio of the 4th order to the 2nd order 

^^""'^ hB JdLo'Nico') 

which at a; = can be approximately estimated for strong classical squeezing (e.g. M^^^ = N) to give IQn^peicPN /hB 
where N is the average photon occupation inside the optical bandwidth. Since this is a perturbative calculation, we 
expect the result, Eq. (P7| to be valid as long as the second term is much smaller than the first one. This will generally 
put a restriction on the intensity, or photon occupation function N(lo), for a given optical bandwidth B. For an optical 
excitation which generates conduction band electrons with kinetic energy of about iOmeV, and an optical energy of 
approximately 1.5eV, we can estimate that for an optical bandwidth KB « lOmeV the average photon occupation 
should be 10"'* — 10~^ for the second term to be 1/10 of the first one {TZ{uj = 0) = 0.1). Usually if one uses correlated 
photons from down-converted light, the spectrum is very wide in the non-degenerate case, and therefore the average 
photon number per mode (N) is very small. For example as discussed in Wang et al.^^ for type-I down converted 
light, a counting rate of lO"' photons/sec was reported over a bandwidth of IQ^^Hz, which corresponds to iV = 10"^. 
In more recent experiments, Bowen et ali^, Heersink et al^^S, and Marquardt et ali^ report of higher intensities for 
polarization-squeezed light. 

It is interesting to note that according to the above estimates there appears to exist an experimentally accessible 
regime where the perturbative approximation breaks down {N > hB/lGir^peid^). In this regime higher orders in 
the light-matter interaction become important, leading to a strong non- linear response to the driving field. It is 
an intriguing direction for future research, especially since strong correlated light sources are becoming increasingly 
availabl o^^'^^ . Another interesting direction to explore is related to effects of non-classical squeezed light \M\ > N on 
the spectrum, which we showed to have a unique effect on the static spin correlations^^. 



VI. EXAMPLES OF PHASE MODULATIONS 



The quantity C{uj) appearing in Eq. is sensitive to the squeezing phases, i.e. the phases oi M{u!), cf. Eq. (|25p . 
In the general case arg{M{uj)) can be expanded in a power series so initially it makes sense to focus on the linear and 
quadratic dependence. First let us rewrite the spectral spin density ([27]) as 

jl^'^ dqz{{S%,^)(^+%^^^^^„ =a{Co^ + C-^fy' [1 + ^Cm (29) 

where Co = oj/B, a — ^^P"^*^ p^' J duj' N{uj'), rj — A-k^ -^^^^^ N , with the electron-hole recombination rate ^rec = 
^ Poptd^ / , the average photon occupation N, and C(a)) — B~^N^'^C{BCj). 

In the linear case we take 9{u} — cjq) — T\ijj — wqI with a phase modulation parameter T = 2 x lO^sec, small average 
photon number N = 10"**, optical bandwidth fiB — 2QmeV, average electron kinetic energy = 30meV, and spin 
relaxation time — 0.2ns. The photon correlation functions were taken to be Gaussian N(oj) — exp[~({uj ^ luq)/ B)'^] 
and M(w) — N(uj)exp[iT\LL! — luqW- One possible experimental realization of a phase modulation of squeezed light is 
by using a pulse shaperiS. With these parameters such a pulse shaper has to turn the squeezing phase a full cycle of 
27T every 1.6GHz. In Fig. [SJa) we draw the r.h.s. of ([29|) divided by a as a function of Cu. We see that the spectrum 
develops as additional oscillatory structure with a frequency T~^, which depends on the type of squeezing (same or 
opposite polarization) . The linear phase modulation in frequency space can be thought of as a time-shift in the optical 
field in time domain. Fig. [SJb). This is somewhat similar to a retardation of part of the random polarization signal, 
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FIG. 5: (a) oj-dependence of the spin structure factor for a linear phase modulation of the squeezing phase of the field (insert: 
the dependence of the phase on the frequency). The figure depicts the 2nd order contribution (dashed blue), and phase sensitive 
4th order contributions for same-polarization (red) and opposite-polarization squeezing (green), respectively, compared to no 
squeezing (dashed black), (b) In the time domain we clearly see that the phase modulation parameter T = 20000s determines 
appearance time of the second correlation peak. The correlation between polarizations determine whether this peak will be 
positive (red) or negative (green). 
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FIG. 6: (a) tj-dependence of the spin structure factor for a quadratic phase modulation of the squeezing phase of the field 
(insert: the dependence of the phase on the frequency). The figure depicts the 2nd order contribution (blue), and phase sensitive 
4th order contributions for same-polarization (red) and opposite-polarization squeezing, respectively, compared to no squeezing 
(green), (b) In time domain, in contrast to the linear case, the chirped phase modulation endows the spin correlation with a 
long tail of correlation determined by T rather than by Ts, and starting immediately after the time Ts. Again the polarization 
correlation determines the sign of the long range temporal correlation. 



leading to the second rise of the spin correlation exactly after that time (T) . This analogy is however not exact since 
the frequency shift is a symmetric function^ ^. 

In the quadratic case (see Fig. [5]), we take 9[uj — o^o) = T{ijj — ujqY with a phase modulation parameter T = 
2 X lO'^sec^, small average photon number N = 10^^, optical bandwidth fiB — 20meV, average electron kinetic 
energy — .iOmeV. The photon correlation functions were taken to be again Gaussian with the squeezing M(lu) = 
N{uj)exp[iT{uj — ujq)'^]- In this case it is clear that the effect of 'chirping' in a;-space, leads through the non-linear 
formula Eq. (j27p to an enhancement (reduction) of the correlations around the DC component of the spin correlations 
when the light is endowed with same (opposite) polarization correlations. In the time domain we see that there are 
two time scales, 77^ and T that govern the behavior of the spin correlations. For opposite polarization correlations 
we observe a qualitative difference: on short time scales t < the correlation is positive but for longer times 
t > T it becomes negative. This change only happens in the opposite polarizations case when the squeezing phase is 
modulated on a scale much smaller than the spin relaxation rate (T~^ ^ 7^). 
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With the experimental parameters estimated above, « 0.1, and a « lO^sec/m^. To estimate the strength of the 
total spin fluctuations we can look at the prefactor a/(a)^ + (js/B)^) for a) = or at the integrated power, depending 
on what kind of experiment is performed. For the first case we get that the prefactor scales like NB'^, while i] cx N/B, 
so the overall strength of the spin fluctuations can be increased while keeping 77 small. Since light with frequencies 
above the e-h gap is strongly absorbed with a typical attenuation coefficient of the order of 10^'^ cm^^ for GaAs, it 
is preferable to assume a thin slab, for example with a thickness of 5/im and a beam of light with an area of Imm? . 
For this geometry and the other parameters given above, the strength of the spin fluctuations around a; = would be 
between 10 — 1000 in units of Ti? ^ with the range reflecting an uncertainty of several parameters. This is the strength 
of the combined contribution of the second and fourth order, so the signal (the fourth order) is estimated from 77 to be 
around 10%. Even though the total spin is very small for a macroscopic sample, we note that recently even a single 
spin has been measured with the Faraday rotation technique^^. 



APPENDIX A: OPTICAL TRANSITION MATRIX ELEMENTS 



We summarize here for convenience the calculation of the dipole interaction matrix elements ^ that appear in 
Eq. ([2])) for optical transitions in direct-gap semiconductors. In the spherical version of the Luttinger-Kohn model'*- 
for degenerate valence bands we can choose the angular momentum quantization axis z to be parallel to the optical 
beam h. This model takes into account the spin-orbit interaction in the valence band, which lead to mixing of the 
heavy-hole (HH) and light-hole (LH) states [j — 3/2). For the fourth order correlation processes which we consider, 
the strongest contributions from LH C transitions are still weaker by a factors of 1/3. For this reason and for 
simplicity we consider only HH ^ C in our model. Let us start with a given electronic wave vector k near the F-point 
in the valence band. The wave function can be written asi 



(Al) 



where fi,m = ±3/2 denote the HH eigenstates, Xm^(^) = T^l^..m\f, 0, 4>) are rotation matrices^^(6', Lp are polar angles 
oik, (j) a global phase), and u"^ are angular momentum eigenfunctions of J^. The quantum number m = ±3/2 denotes 
the angular momentum projection in the direction of k, and fi = ±3/2 is the projection in the direction z. Using the 
usual notation for the wave functions of s and p orbitals, Ug^j — ^ ^ ^^'^ ^-3/2 ^ ^ ^" ^'-'^ 

conduction band V'fe s ~ e^'^^Ug, where u'^^^ = S" T "^^-1/2 ^ ^ i- The electronic part of the dipole matrix elements 
are then given by 



(A2) 



where I? = ef* is the dipole moment operator. Let us consider the matrix elements for the transition into the state 
fc,s = 1/2 



^k,k' 



Sk.k' 



-3/.,+3/.(fc)(^ ^\^\-^(^ + ''^) + X-3/.,-3/. 



V2' 



(k){S } \D\^{X - iY) i) 



-D 



= Sk,k'X + 3/2, + 3/2ik)—/^ix + 



4,fc'X_3/2,+3/2(fc) — + 



T2' 



where we used the facti that the only non-zero matrix elements are {S\Dx\X) = {S\Dy\Y) = {S\Dz\Z) — D due 
to the spatial symmetry of these integrals. We see that due to spin orthogonality only one of the orbital angular 
component, fjL = +3/2 contributes to the transition into +1/2 state (a similar expression can be derived for s = —1/2). 
The full interaction matrix elements are given by d'^'l'i"!^; = £ujf._^,^\ ■ (-D)^'™, where A = ± denotes the helicity and 
^i^k^k' ^ \/^^k-k' /'2'£qV. We see that A is determined by the conduction band spin, i.e. s <4> A, while for an arbitrary 
direction of k both m = ±3/2 contribute to the transition. We can therefore omit the index A and summation over 
it, and continue with a simplified interaction Hamiltonian 

H d-k^'(i^p,sbp,s'Vk,m + h.c. (A4) 

k,p s,m 
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FIG. 7; (Color online) A schematic for the time evolution for (a) direct evaluation of the static correlations via second order 
perturbation theory (b) the two-times spin fluctuation function with a dissipative part shown in grey, and (c) the evaluation of 
the same correlation with equations of motion method. 

where we now use s also for the helicity. By applying a transformation to the basis where angular momentum 
projection is along n with new operators Vk.fi 

Vk,m = X„,+3/2 (^)Wfe,+3/2 + X„,_3/2 ik)vk,~3/2 (A5) 

the interaction appears as 

fc,p s,m 

= J2J2'^1+P.sh'' [K;r^^^ + 3/2, + 3/2(fc)^^fe,+3/2 +<;p^^\ + 3/2,-3/2(^)^'fc,-3/2) + 

k.p s 



+ «;p^^^X-3/2.+3/2(fc)^^fe,+3/2 + <;p^^^X-3/2,-3/2(^)^'fe,-3/2) 



h. 



c. 



^^""l+pjp^^ [K;r^^X + 3/2.+3/2(fc) +<;p^^^X_3/2, + 3/2(fc))^^fc,+3/2- 

k.p s 



+ (c^fe,r^^X + 3/2,-3/2(^) + t^fe,p^^^X-3/2,-3/2(fc))^'fe,-3/2 



h. 



c. 



The rotation matrix element can be written aa^ 2?^, mi^^ — '^'^m' ml^)^™"^ (here d is not the dipole matrix 
element). By using the symmetry relation d^^, ^{9) — (—1)™ -mi^) ^^'^ setting = 0, it is straightforward 

to show that the coefficient in front of Vk,+3/2 vanishes when s — —1/2 and the coefficient of Vk^-^3/2 vanishes when 
s — +1/2. Therefore there is a one-to-one correspondence between the projection of angular momentum of the initial 
valence band state 11 and the projection of the spin of the final conduction band state s. Renaming v as v, we can 
rewrite the interaction term as 

J2T.'i+P,s^'^^-')^^-MlM+xl-sik))vk,. + h.c. (A7) 

k,p s 

where a common index s = ±1 <^ ±1/2 <^ ±3/2 is used for angular momentum, spin, and helicity. In cases where 
the angular dependence does not lead to other consequences other than some numerical prefactor (in cases where an 
isotropic approximation for all the other parts in the integrand is reasonable), we just omitted it and assume that D 
absorbs the prefactor. 

APPENDIX B: PHENOMENOLOGICAL LEVEL WIDTHS 

The fourth order contribution {S{t)S{t'))'}^l which we calculate here should match our previously calculated in- 
jection rate of static correlations^^ in the limit t' t. Indeed it can be shown that they give the same result if in 
the previous calculation a finite level width is incorporated for the finite state propagator, and the intermediate 
state level width is taken to zero. To explain the source of this problem we need to compare the two approaches. 
Previously we calculated the static correlation (t = t') directly, using a second order perturbation expansion of the 
wave function. In interaction picture the static correlation can be written as (see Fig. [7fa)) 

( {^el\U\t, -00)5(t)2C/(i, -00)|^Aei) ) f^eld (Bl) 
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while here we calculate the two-times correlation, which is defined as (see Fig. [7jb)) 

( {4'el\UHt,-^)S{t)oUit,t')Sit')oU{t',~CX^)\^Pel) )freld (B2) 

where S{t)o denote operators in the interaction picture. Since a direct evaluation of (|B2[) is very difficult for fourth 
order processes, we use the equations of motions for S{t), S{t'), which means actually that we split the intermediate 
time evolution into a product which gives (see Fig. Ulcj) 



—oo) X 

U\t', -CX^)S{t')oU{t', ) f^eld 



(B3) 



which is in principle equivalent to (|B2j) due to unitarity. The difficult issue is how to account also for the relaxation 
processes in such calculation in the simplest way while retaining consistency. In our calculations we added them 
explicitly as level broadenings (7). However when the relaxation is added to the time evolution it is no longer possible 
to claim that unitarity applies 



(B4) 



and therefore in principle different approaches can give different results. Specifically for the limit t t' the evolution 
depicted in Fig. [TKa) might not give the same result as going along [t — > — cxd — > t) in Fig. [Tfc). In principle this shows 
that adding dissipations phenomenologically in higher order perturbation theory is very tricky and it would be better 
to have a microscopic model instead. We have shown^ another approach with Keldysh diagrammatic formalism that 
is promising to lead in the future to a derivation in which the dissipation processes will be derived from a microscopic 
Hamiltonian, and will not suffer from such ambiguities. 



APPENDIX C: CALCULATION OF THE {S'^^\q,t)Sf\q' ,t')) CONTRIBUTION 

We begin by writing explicitly the second contribution to the spin wave operator 

5,(<?,t)(2) =-|d|2e-^=* f dhe-^'^e-^^'^ f'd^e^'- ^ [ (CI) 

■^*o ■^*o k,p,p' 

e 9. p e ^ " f p ' Op,| (nfc_^ fc„p+p, |(t2jo - n-fe_g_p',fe_p,|(,i2joj + 

_ i(Acfc_,_fc_p-ajp)ti -i(Aefc_,,fc_p-cj ,)t2L* ( c (fn^r. — (fn"),^) h ,— 

_p-j(Aefc,fc_p-5-Wp)tl i(Aefc,fc_p_,-w /)t2L* / c Jfn)n—nV , Jfn\n) h n 

^ CP >i V'-fc-p-g+p',fe,il''2j0 "-fe-p-g,fc-pM l>'-2j0j Op'i 

In order to derive the third order contribution Sz{q, t)'-^-' it is necessary to formally solve the equations of motion for 
the various operators that appear in the second order contribution (jCip . We start by writing integral expressions 
for n^'^, , again employing the basic equations of motion ([7]), and repeat the approximations of static spin waves, by 
neglecting the free evolution of Uk^k' terms 

''^k—q,k—p+p',a{^2) 

= e-^*^ / ' dtse-^'^Y, [-rfnLp+p,_p»,,_,,.(i3)oV'.fe)o + d*b;„At3)oPk-q-p",k-p+p'At3)o 

■'to p„ 



^e-^*^ dtse^'H [-e^(^^-— +-'-"-"''"^*■'rfP^p+p,_p.,,_,,,V'.- 



+e-^i^->^-P+p',>^-.^p"-^p")'^d*b;„^Pk^,.p„,k.p+p,,^ 



E 

p" 



k — q.k — p-\-p —p" 



k—p+p'—p",k—q,a P 



7 + i{Aek~q,k-p+p'-p" — ujp") 

p"cr- 'k-q-p" ,k-p+p' ,tT 



k — p-\-p,k — q — p 



7 — i{Aek-p+p' ^k-q-p" — Wp") 



(C2) 
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'^fc-p+p'-g,fe,<7(*2) — 

J to p// 

._// -J tn 



E 



_e»(^«fe-P+p'-9,<=-p"-'^p")t2(;pt 5 

k—p",k—p+p' — q,a P " 

7 + i{Aek-p+p' -q,k-p" — ijJp") 



+ 



+ 



e-i(^^>^.^-.+.'-p"-,-'-p")t-d*b;„^Pk-p+p>-pn-g,k,a 

7 — «(Aefc,fc_p+p'_p"_g — LOp") 



(C3) 



'^fe-p'-g,fc-p,cr(*2) 



-7*2 



/ (^ise''*^^ E [f^'Pfc-p,fe-p'-g+p",<7(*3)oVv(i3)o - d*b;„^{t3)oPk-p'-q,k-p+p",a{t3)o 

J to pll 

-e-i(^'>>-P+p".>^-p'-.-^."^'^d%;„,Pk-p'-q,k-p+p",a 
oii^^k-o> -a+n" .k-t,-<^n") 

E 



f.i(^<^k-p'-,+p",k-p-<^p")t^CiPl u , , „ K"a 

k—p,k—p'—q+p",(7 P 



7 + i{l^ek-p' -q+p" ,k-p — <^p") 



e-i(^^>'-P+p",'^-p'-,--'p")^^d*hl„„Pk-p'-q,k-p+p" 
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(C5) 



When now substituting n^'li in the expression for Sz {q, t), we only need the oc terms, because in the correlation 
1x3 only averages of the form (PP^) are contributing for zero temperature, leading to 



•'*<' k,p,p',p" 



(C6) 
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FIG. 8: Diagrammatic representation of the 1x3 excitation process with correlated photons. The diagram shows the emitted 
spin- wave (wiggly lines to the left & right of the diagram), pump light correlations (wiggly lines inside) and particle- hole 
excitations (solid & dashed lines, respectively). 
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The contraction {S^\q,t)S^^\q,t)) leads to contractions such as 



(6*6*,6p"6p) (Pfe_p_^ feP^-p+p'-p" 



(C7) 



which result in the constraint k — k — q which restrict the phase space of fc,fc summation considerably, compared 
to the 2x2 contribution. Diagrammatically, these fourth order contributions can be described as one fermion loop 
diagrams, see Fig. |H1 
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